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ABSTRACT 


We consider the constraints that the available X-ray spectral and imaging 
data place on the mass distribution and mass to light ratio of rich 
clusters. We find for the best determined cases that the mass to light ratio 
is less than 125 hg 0 at radii exceeding 1 hgQ Mpc. The mass to light ratio is 
approximately constant at radii exceeding 1 h 50 Mpc but may rise to values of 
roughly 200 hg 0 in the central regions. The fraction of the total mass that 
is in baryons, primarily the hot X-ray emitting gas, is roughly 3035 thus 
setting the mass to light ratio of the "dark" material to roughly 70. The 
model that fits the X-ray data for Coma is in good agreement with the observed 
optical velocity dispersion vs. radius data. 


Subject headings : cosmology - galaxies: clustering - galaxies: X-rays 
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1. Introduction 


It has become a part of the astrophysical lore to suppose that the 
virial masses in the rich clusters are large. Typical mass to light ratios 
(M/Ly, where Ly Is the visual luminosity) are generally supposed to lie 
between 200 and 300 h'^Q* (for a recent review see Rood 1981). This missing 
mass (or missing light) problem fits neatly into a strongly held picture that 
the mass-to-light ratio Increases with scale size of the bound system, and 
that the universe Is closed, as required by the Inflationary cosmologies. Of 
course, even the large mass to light ratios of clusters of 300 h 50 Is too low 
a value to provide the closure M/L ratio of 1000 h 50 (Felten 1985,1986). 
However, as is well known, the basis for these views are very fragile. 

Detailed optical data on light profiles and velocity dispersions as functions 
of radius are available for only a few clusters (c.f. Kent and Gunn 1982, Kent 
and Sargent 1983). Even where such data does exist translation into a virial 
mass profile requires many assumptions, including spherical symmetry, constant 
mass to light ratios, details of the distribution functions of the galaxies 
(the form of the orbits e.g. radial vs. symmetric, the assumption of 
relaxation of the galaxies in the potential) and the form of the potential. 

The final answer is quite sensitive to many of these assumptions. For 
example, Bailey (1982) has shown that breaking the constant M/Ly requirement 
and allowing the virial mass to be centrally concentrated allows the total 
mass to be reduced: M/Ly could be as low as 50 h~* 50 in the Coma cluster in 

such models. If we can demonstrate that this type of model Is preferred, the 
impact on our cosmological viewpoint should be profound. 


*h 50 is the Hubble constant in units of 50 km s - * Mpc - * 
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One of the ways to test this question is through study of the X-ray 
atmospheres of the clusters. Determining total mass to light ratios in 
clusters using observations of the X-ray emitting gas should hold many 
advantages over traditional optical methods. The gas is guaranteed to have an 
isotropic velocity distribution and to be quasi-static in the inner cluster 
(Cavaliere 1980). Thus integration of the hydrostatic equation can allow a 
direct measurement of the run of the total cluster mass. The reason such 
determinations have not been seriously pursued, as yet, is that it is 
generally assumed that both gas density and temperature profiles must be well 
known to use the method (e.g. Sarazin 1985), and at present only the gas 
density profiles are available from the Einstein Imaging data (Jones and 
Forman 1984). There have been attempts to measure the mass of galaxies by 
this method most notably that of M87 (Fabricant and Gorenstein 1983). 

However, the density profile data, when combined with integral spectral 
information, can provide very strong constraints on cluster virial masses even 
without the full temperature profile. 

The emission measure weighted spectra for many clusters are available 
from 0S0-8, Ariel-V and HEA0-1 data and for the past few years these spectra 
have been approximated by isothermal bremsstrahlung emission (e.g. Mushotzky 
1984). However, recent analyses (Henriksen and Mushotzky 1986b) have shown 
that such single-temperature models are not good fits to a number of 
clusters. In these cases a range of temperatures is required, including 

material at both higher and lower temperatures than would be inferred from an 
isothermal fit. Henriksen and Mushotzky have analyzed this data in terms of a 
polytropic equation of state which allows a simple two parameter fit to the 
data. They then used this fit to measure virial masses in clusters with the 
hydrostatic equation, and concluded that the virial masses were considerably 
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lower than previously estimated. 

After summarizing the necessary background material in section 2, we 
generalize this procedure in section 3 and show that within such a polytropic 
model a value of the mass to light ratio M/L v ~ 125 h 50 - ^ at R = 2 h 50 -1 Mpc 
for both Coma and Perseus clusters is entirely consistent with all the 
available optical and X-ray data. 

However, the assumption of the polytropic equation of state may be 
considerably too specific, and in section 4 we consider an alternative 
procedure. In this section we assume that the mass-to-light ratio is constant 
in the clusters (an assumption familiar from previous analysis but still not a 
necessary assumption) and that the light profile is given by an analytic King 
model (which except in the central regions of Coma seems to be an adequate 
representation of the data (Quintana 1979)). Knowing the virial mass profile 
and the density profile, we can use the hydrostatic equation to obtain the 
temperature profile as a function of three parameters - the mass to light 
ratio, the optical core radius, and the asymptotic temperature. We can then 
construct the integrated spectrum. We find that this model cannot provide an 
adequate fit to the spectra, but invariably gives too isothermal a temperature 
profile. Inclusion of a large central mass in the cluster core (~ 2 x 10 13 
M Q ) allows us to obtain agreement with the spectral data, but when this is 
done the required M/L ratio outside the core again falls to low values (< 125 
h 5 g "1). Preferred values are typically around 60 h^Q - *. Finally, we show in 
section 5 that there are models which provide an adequate fit to the optical 
as well as the X-ray data, and that these models generally give M/L v » 100 
h 50”*» be y° nd 0*5 h 50 Mpc. 

In section 6 we summarize our conclusions, and consider the implications 


of our results. 
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2. Background Material 


a. Gas Density Profiles 


The Imaging Proportional Counter (IPC) aboard the Einstein Observatory 
obtained X-ray surface brightness profiles for many clusters out to radii of 
roughly 1 Mpc (c.f. Abramopoulos and Ku 1983, Jones and Forman 1984). For a 
spherically symmetric atmosphere the surface brightness S(r) at projected 
radius r is related to the emission per unit volume into the IPC passband E(R) 
at radius R by the equation, 


S(r) 



Ef R) RdR 

(R 2 - r 2 ) X ' 2 


This Abel equation can be inverted to give 


(1) 


E(r) 



_d 

dr 


T R S (R) dR 

Jr (R 2 -r 2 ) 


( 2 ) 


as noted by Cavaliere (1980). Direct inversion of equation (2) without 
additional external constraints is quite unstable, (the deconvolution 
technique used by Fabian et §2. (1981) and Stewart et aj. (1984) is equivalent 
to inversion of this equation with additional constraints on the potential, 
total X-ray flux and mean X-ray temperature). The simplest procedure, fitting 
a smooth function to the data prior to inversion, does retain most of the 
available information in the surface brightness profile. As Jones and Forman 
have shown, the surface-brightness profiles are well fit, outside the very 
central core, by a law of the form 
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2 

S(r) = S n (1 + ) 1/2 - 36 (3) 

a 

This functional form was originally derived for the case of an isothermal gas 
and isothermal cluster (Cavaliere and Fusco-Femiano 1976), but this 
interpretation is not consistent with the optical velocity dispersions, 

optical size properties and the X-ray spectral data (Mushotzky 1984). As will 
be clear from our later discussion, equation (1) is best considered as an 
extremely good empirical fit with the surface brightness described by equation 
(3). Equation (1) may now be inverted, using equation (2) subject to the 

assumption of spherical symmetry, to give an emission per unit volume 

E < r > = 4 - 1/2 W-- 1/2) < 1 + ir>' 3s < 4 > 

E(r) * n e 2 e (T) where e(T) is the emissivity of the gas convolved through 
the Einstein IPC passband (and r is the gamma function). Numerical 

integration of a thin bremsstrahlung spectrum through the effective area of 
the IPC as a function of energy shows that e is a weak function of temperature 
(a t~ 0*35) (e.g. Fabricant, Lecar, and Gorenstein 1980; Fabian et al. 1981) at 
temperatures greater than 2 keV. Thus e is almost entirely dependent on 
density. In particular, if the gas is described by a polytropic equation of 
state T ~ n Y “l, then the density profile may be written 


n - n Q (1 + r 2 /a 2 ) 


(5) 


where 


6 = ~ y ~ 6 / (1 - 0.18 ( Y - 1 )) 


( 6 ) 
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Jones and Forman found a maximal range of 6 of .40 to .83 (with a minimal 
range of .52 to .68) in their clusters which translates to s = 0.60 1.25 for 

y = 1 and 6 = 0.54 -*• 1.11 for y = 5/3. The 6 ranges and other X-ray 
parameters for the five clusters we shall consider here are summarized in 
Table 1. 

TABLE 1: Polytropic Model Parameters 


Cluster 

6 

T IS0 

T o 

Y 

A85 

0.9-0.98 3 

6.8±0.5 C 

12.8-22. 5 e 

1.3-1.66 

A426 

0.83-0.90 a 

6.4±0.4 C 

7.6-10.8 e 

1.15-1.36 

A1656 

1.00-1.30 b 

8.0±0.4 d 

15.3-20 d 

1.38-1.6 

A1795 

0.98-1.20 3 

e + 0 . 5 C 
6,5 -0.3 

10.2-21 

1.25-1.66 

A2199 

0.95-1. 10 e 

3 . 6±0 . 6 C 

4-9.5 

1.0-1. 5 

where T Q 

is the central 

temperature 

in keV 



Tj S q is the temperature if the gas is isothermal 


a Jones and Forman 1984 
b Abramopoulos, Chanan and Ku 1981 
c Mushotzky 1984 
d Henriksen and Mushotzky 1986 
e Henriksen and Mushotzky 1986b 
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As Forman and Jones note a number of clusters are not well fitted by 
equation (4) within the core regions. These are generally the clusters which 
have cooling flows inside the core. Clusters, such as Coma, which do not have 
such flows are generally well fitted throughout. It is possible, that the 
best fit value of y may be changed by the presence of the strong cooling flow 
in A85, A426, A1795 and A2199. However, it is unlikely that the presence of 
the cooling flow will change the value of the derived central temperature. 
Given that the cluster without a cooling flow. Coma, has the highest allowed 
value of y, it may be possible that the cooling flow to some extent hides the 
signature of a higher y polytropic model. 

The density deconvolution outside the core does not, of course, depend in 
detail on the core properties. The density deconvolution is also not 
sensitive to the details of the fit to the surface brightness or to the 
symmetry assumption. In particular, as Fabian et al. (1981) and Rybicki, 
Gorenstein and Fabricant (1984) have shown, if the clusters are cylindrical ly 
symmetric and the ellipticity is not too high the density profiles do not 
differ radically from those obtained assuming spherical symmetry. Further, 
much of this weak symmetry dependence cancels out in virial mass 
determinations. 


b. Temperature Information 

The integrated emission from the cluster is dominated by emission from 
the regions where the density profile is determined by the Einstein 
measurements. Temperatures based on the isothermal models for the HEAO-1 data 
are given in Table 1, but such models do not generally provide an adequate 
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fit. More generally, the spectra may be fully specified by the differential 
volume emission-measure as a function of temperature. A function of the form 


E(T) dT = A {-!-> U [1 - { 4 -> W ] 

'o 'o 


1/2 


dT 


(7) 


does provide an acceptable fit to the HEAO 1 data. Equation (4) is obtained 
for the particular case of a polytropic equation of state — I — = (— — ) Y “ * 

'o n 0 

and the density profile of equation (5) with 


A = 


? J n 2 

Zira riQ 


(8a) 


u = 


(3-y) 6 - 1.5 


(8b) 


W <f> 


(8c) 


where <t> = «(y - 1). The value of w is extremely poorly constrained by the 
observations; if it is force fitted to a typical value of ~ 3 then the values 
of T Q and y(where the ranges indicate 95% confidence limits) for a number of 
clusters are those given in Table 1. T Q is essentially a monotonically 
increasing function of y • Values of T Q are larger for the case 6=1. 

The preferred temperature profiles are generally quite shallow and in 
some cases (e.g., A2199) could be close to isothermal. Coma has one of the 
steepest slopes and here the temperature must vary by about a factor of 2 - 3 
within the radius where 90% of the flux originates. For the polytropic case 
with y = 1.4 and 6 = 1, the temperature in Coma falls from 17 keV at the 
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center to 8.8 keV at 2 core radii and 4.6 keV at 5 core radii, if y = 1.25 the 
temperature varies from 15.3 keV at the center to 10.2 keV at 2 core radii and 
6.8 keV at 5 core radii. An isothermal model by contrast (Table 1) gives a 
best-fit temperature of ~ 8 keV (Henriksen and Mushotzky 1985). 

c. Hydrostatic Balance 


The hydrostatic equation in a spherically symmetric cluster may be 
written as (e.g. Fabricant, Lecar, and Gorenstein 1981) 


M v (r) 


kTr 
G m 


r d In n d In T i 

1 d In r d Inr 1 


(9) 


where My is the virial binding mass as a function of radius, T is the 
temperature, n the density profile and i is the mass per gas particle. With 
a density profile of the form of equation (5) this may be rewritten as. 


M v (r) 


ka 

Gm 


Tx 


2 

26 

(1 + X 2 ) 


+ 


d In T , 
d In r J * 


( 10 ) 


where for simplicity we have defined x = r/a, where a is the gas core 
radius. For a polytropic equation of state, this simplifies still further to 


M 


V 


2 Y 6 


ka T o 
G m 


(l+x 2 ) 


!+<)> 


(11a) 


where T Q is the central temperature. Putting in typical values for the 
parameters and substituting we find 


M v = 1.92 x 10 14 6 Y (T Q /10) keV (a/0.25) Mpc x 3 (1 + x 2 )" (1+ * } M Q (llb) 
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The ratio of the "polytropic" to "isothermal" binding mass is 


Vl> - y x -2* 


(lie) 


which for "typical" values of the parameters (see section (3)) is 


M v (y) 

M v (iso) 


- 1.25 x 


- 1/2 


We shall turn to this simple case first in section (3) before proceeding 
to a more general description in section (4). 
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3. Polytropic Equation of State 


a. General Description 


We must always recognize that a polytropic equation of state corresponds 
to a particular choice of temperature profile which may force us to 
conclusions which can be avoided in a freer model. However, the 
simplification of such a parameterized fit allows analytic descriptions and 
gives considerable insight. 

Equations (7) and (11) provide the basic equations for this section. For 
a given cluster we may use the following procedure. First, the range of 
acceptable 6 is determined from the IPC data on the cluster. Generally 6 lies 
around 1 and fits extend from slightly less than 1 h 50 1 Mpc to about 2 1/2 
h 50~* M P C ( Jones and Forman 1984). For each value of & a range of y and T Q 
( Y ) is determined, using equation (7) in conjunction with the HEAO 1 data, as 
summarized in Table 1. Finally the run of virial mass and the range of the 
determinations can be found from equation (11). 

Before proceeding to analyze specific cases, some general points follow 
at once from the functional form of equation (11). The density corresponding 
to this mass profile is given by 


p 


y6 


2irGa 



3 - [2 » - 1] 

(1 + x 2 )<* + 2 ) 



( 12 ) 


From equation (12) we can see that for 2$ i 1, the virial density falls to 
zero at a finite value 


x max 


3 ,1/2 

2 <> - 1 J 


[ 


(13) 
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A cutoff in the virial density at such small radii (typically about 2 core 
radii) is not acceptable since both gas and galaxy light extend beyond this 
point. Thus, we must impose the condition 


* s 1/2 


(14) 


so that there is no "natural cutoff " at a few core radii. 

This constraint should be viewed as empirical relation. Gas obeying the 
form of equation (5) must have a low polytropic index to extend to even 
moderate radii. 

As can be seen from Table 1, the physical constraint on <t> favors lower 
values of 6 and shallower slopes (smaller y). Coma pushes hardest on this 
constraint. The best fit 6 (1.15) and lowest allowed value of y (1.4) would 
give 4> = .46 which would result in a very rapid mass drop beyond the core 

radius. 

If we adopt an analytic "King" approximation to the cluster light profile 

»L ' » Lo d + -jj-) ' 3/2 • (15) 


(where is the density of light from the cluster) then we may write the M/L 
ratio as a function of position 



I ll ll 

2nGa^ p^ Q 


kT 


m 


3 - ( ? *-l)x 2 , (1 + B 2 X 2 ) 3/2 


(l+x‘) 


2 * 6+2 


(16) 


where B 2 = a 2 /b 2 


is the square of the ratio of X-ray to optical core radii. 


As opposed to previous treatments we do not require the mass profile to follow 
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the "King" light profile and we adopt this form for ease of comparison with 
previous work. 

b. Global Mass and M/L 

The virial mass at larger radii is maximized if we adopt the lowest 
acceptable values of 6 and y and the largest acceptable X-ray core radius 
(Table 2). The maximum mass profiles for the five clusters of Table 2 are 
shown in Fig. 1 together with the corresponding M/L v profiles for Coma and 
Perseus. In computing the mass-to-light profiles we have used King model fits 
to the light distribution of the form of equation (15) with the optical core 
radius b = 0.37 h"* 5 Q M p c and Ly inside 7 hgg"* Mpc = 1.6. 10 + * 3 h 50 “ 2 L Q for 
Coma (Kent and Gunn 1982) as modified in Kent and Sargent (1983), and b = 0.34 
h 50~* ^pc and L V ins ^ de 5.6 h 5 Q-l Mpc = 1.1.10*3 h 50 -2 L Q for Perseus (Kent 
and Sargent 1983). The corresponding visual luminosity profiles are 


L v = 

6.2.10 12 h 5 „- 2 f(y) 

(Coma) 

(17a) 

L v = 

4.4. 10 12 hg 0 -2 f (y) 

(Perseus) 

(17b) 


where y = r/b and 

f (y) = In (y + (1+y 2 ) 1/2 ) - y/ (1+y 2 ) 1/2 (18) 

The visual luminosity inside 2h 50 “* Mpc is 8.7 10* 2 h 50 -2 L Q (Coma) and 6.5 
10* 2 h 50 -2 L Q (Perseus). Abell (1977) gives the visual luminosity in Coma out 
to a radius of 8.6 hgg - * Mpc as (1.6-2.75) 10*3 |_^ Using his normalization 
would change the normalization of the total luminosity in Coma to (6.8-11.7) x 
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10 12 rather than 6.2 x 10* 2 L Q and result in even lower values of M/L v . 

TABLE 2: MAXIMAL BINDING MASSES 


Cluster 

Core Radius Range 3 

Maximal Mass 

A85 

0.19 - 0.26 

6.7 

A426 

0.23 - 0.34 

7.3 

A1656 

0.36 - 0.47 c 

9.3 

A1795 

0.20 - 0.40 

8.5 

A2199 

0.12 - 0.16 

6.7 


a in Mpc from Jones and Forman 1984 
^at 2 Mpc in units of 10^ M Q 
c Abramopoulos, Chauan and Ku 1981 

Even with these extreme parameters (lowest allowable y and 6), the M/L 
ratio (using Kent and Gunn's values) is only approximately 125 h 50 - ^ for Coma 
and Perseus within 2 h 50 "^ Mpc compared to Kent and Sargent's (1983) best 
guess values from the optical data of 160 h 50 ~* for Coma and 300 h 50 -1 for 
Perseus. With a more probable choice of parameters see section 5 (e.g. y = 
1.37, T Q = 15 keV, 6 = 1.14 for Coma) the values are considerably lower, with 
best fit values of ~ 90 for M/L y in Coma at 2 h 50 -1 Mpc. An interesting 
point about the mass to light ratios is their peak within the core followed by 
a relatively constant value at larger radii. We shall show in the next 
section that this is not simply an artifact of the polytropic assumption but a 
fairly direct consequence of the X-ray data combined with the assumption of a 
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King light profile in the central regions. It is possible that the innermost 
regions of clusters are not well described by a King model (Quintana 1979, 
Beers and Tonry 1986). To see whether the rise in M/L in the central regions 
of Coma is due to the artifact of an analytic King model assumption for the 
light distribution we take the result of Quintana who shows that while in the 
central O^hgQ -1 Mpc there is an excess of light of 1.5-2 times the King model 
density, the King model is a good fit at larger radii. This effect would be 
to put an inflection in the mass to light ratio curve with the central value 
of M/L dropping to ~ 120 (see Fig. 4), the value seen at radii greater than 1 
hgQ -1 Mpc. ^ we 1ns ^ s ^ that the M/L curve be smooth, with no inflections in 
M/L vs. R, then we can speculate that the true M/L vs R curve resembles the 
dashed curve in Figure 4. It is thus entirely possible that there is not a 
rise in M/L in the central regions of clusters if the galaxy distribution is 
steeper then a King profile. 


c. Total Mass In Baryons 

The mass in gas can be analytically calculated for 6 = 1 as 
Mga $ (< x) = 4.5xl0 13 M q (a/0.4 Mpc) (n Q /2xl0' 3 ) (x - tan' 1 x) h 5Q ~ 5/2 (19) 

O 

where n Q is the central gas density in parti cl es/cm. For Coma at 2 Mpc 
Mg a$ - 1.6 x 10 14 M q . (If we use y = 1.15 Mg as = 1.16 x 10 14 M Q ). If we use 
for the mass in stars a M/L v ~ 7, appropriate for an elliptical galaxy 
(Pickles 1985), then the mass in stars inside 2 h 50 _1 Mpc is M* - 6 x 10 13 M Q 
giving a total baryonic mass of 2.25 x 10 14 M Q at 2 Mpc. The baryons thus 
contribute at least 20% (16% if y = 1.15) of the total mass inside this 



18 


radius. Therefore the M/L y of the "dark matter" in Coma is < 99. As an 
extremum for the "minimum" allowed virial mass and M/L y we take the <t> = 0.46 
case which gives a total M/L y for Coma of 44, a mass of 7.2 x 10 14 M Q for 
Abells largest allowed luminosity, and a M/L v for the "dark matter" of only 
~ 30. It is thus not inconsistent to state that the dark matter in clusters 
(at least Coma) maybe quite similar in form to that postulated for the disk of 
spiral galaxies (van Albada, Bahcall, Bergeman and Sancisi 1985). 

4. Generalized Models 


An alternative approach to the problem is to assume a mass profile and 
then to use the hydrostatic equation to determine the temperature profile 
(this method is similar to that of Sarazin and Bahcall 1977). One may then 
test if the integrated spectrum is an acceptable fit to the data. This method 
has the "virtue" of making no a priori assumptions about the functional form 
of the temperature profile. 

For this case it is simplest to use the hydrostatic equation in the form. 


T = (1 + x 2 ) 6 


Tr + mG 
(1+x 2 r ) 6 ka 


x R M v (x) 

Jx x 2 (l+x 2 ) 6 


dx 


( 20 ) 


where T R is the temperature at a reference radius x R . 

First, consider the case of constant M/L ratio models where the 

luminosity has the analytic King model form of equation (18). This 
integration is shown in Fig 2 for Coma with a = 0.34 h 50 ~* Mpc and b = 0.37 
h 50~* Mpc ^ or cases 6 = and $ * 1.25, for M/Ly = 100 and 150, and for 

T r = 0, 10 and 20 keV at 2 hgQ - -'- Mpc. 

The integrated spectra of these models appear very nearly isothermal 
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since the temperature profiles in the core regions are much shallower than 
those of the best fitting polytropes. If the integrated spectra are force- 
fitted to the HEAO-1 data, they give M/Ly * 110 h 50 -1 for 6 = 1.25, but they 
are not acceptable fits. 

In order to obtain a better fit to the spectrum we need a steeper 
temperature profile in the central regions, which in turn requires additional 
mass in the core regions. Inclusion of such additional mass produces mass 
distributions similar to those of the polytropic models of fig 1, with high 
mass-to-light ratios in the core and relatively constant values outside. The 

1 O 1 

required excess core mass in Coma is about l-2.10 Aa h 50 “ A M Q and the mass to 
light ratio is about 100 h 5 Q _ * beyond the core. 

5. Optical Velocity Dispersions 


So far we have not considered whether the models of the previous two 
sections are consistent with the optical velocity dispersion data. It is 
simplest for this purpose to use the analytic form of equation 11 which we now 
write in the form 


M, 


M_ 


(1 2 + x) (1+ * } 


( 21 ) 


where M = 26y - — - and x = r/a. Now assuming the velocity distribution 
0 Gm 

of the galaxies is isotropic, and the galaxies are in hydrostatic equilibrium 
in the potential well of the cluster, the light weighted projected velocity 
dispersion at projected radius p (where p is in core radii) is given by 
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(P) - 


I, 


(1+x 2 ) 1+ * 


(x 2 - p 2 ) !/ 2 
fl+B 2 x 2 ) 3 ' 2 


dx. 


j: 


x d x 


P (x 2 - p 2 ) W (1+B 2 X 2 ) 3 ' 2 


which for Coma, where B = 1, may be simply integrated to give 


/ IT 


(2* + 3) rlt+3/4) 


Y 6kT 


o /i -<t> 


(i+r) 


(22a) 


m 


The ratio of the r functions is approximately constant over the interesting 
domain of <j> and we can thus write 


a p (0) * 866 ( Y 6) 1/2 (T/10)J/ 2 km/sec (22b) 

As can be seen from fig 3a, the models at the low end of the Y and 6 values 
determined from the X-ray data provide a very good fit to Coma's velocity 
dispersion profile while higher values do not. The allowed values from 
optical data only for Coma have 6.8 < T q < 15.8 and 1.1 < Y < 1.34 and thus 
only marginally overlap the X-ray best fits (Henriksen and Mushotzky 1986). 
"Isothermal models" do not provide a good fit either. Of course the 
appropriate "King" models (Kent and Gunn 1982) can fit the velocity data. The 
intersection of the best fit values for both optical and X-rays is small and 
rather marginally acceptable as regards the X-ray data. If we fit the Coma 
optical data within the region where the x-ray data provide any constraints (< 
9 core radii or 3.6 Mpc) then the optical data can be well fit (reduced chi- 
square < 1.5) with parameters (6 = 0.9, Y = 1.4 and T(0)= 18 keV) which are 
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entirely consistent with the x-ray data alone. Using the largest y and T Q (at 
68# confidence) from the optical data alone we obtain the mass-to-light 
profile shown in Fig. 4. The mass-to-light ratios have the extremely 
interesting property of being constant at larger radii and have asymptotic 
values of 100-125 h 5 g _ l. 

None of the models provide a particularly good fit to the data in Perseus 
- the central velocity dispersion is always too high (Fig. 3b). This is a 
notorious problem in Perseus and is usually attributed to anisotropy of the 
galaxy velocity distribution in the core (e.g., Kent and Sargent 1983 and 
references therein). Because we have assumed isotropic galaxy orbits the 
magnitude of the temperature discrepancy is almost exactly the same as that 
found by Kent and Sargent (1983), a factor of two. We note that if we 
renormalize the central velocity dispersion by this factor the predicted trend 
of velocity dispersion with radius follows the optical data quite well. 
Another explanation might be the presence of a foreground clump of high 
velocity galaxies projected on the core of Perseus. 

The only other cluster which has sufficiently high quality optical data to 
attempt a fit even for o p (o) is Abell 2199. In this case the combined optical 
and X-ray data have small allowed boundaries (Henriksen and Mushotzky 1986b) 

with 5.0 < T < 5.8 keV and 1.1 < y < 1.2. 

~ 0 ~ ~ 

The general agreement of optical and X-ray data for Coma and Abell 2199 
seem to rule out gross anisotropies in the orbits of the galaxies which, 
however, seems to be required for Perseus. 
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6. Conclusions 


Our final conclusion is that the X-ray data cannot accept high mass to 
light ratios in these clusters. Typical mass to light ratios must be less 
than 125 h 50 '* on the basis of the X-ray data and preferred values are 
considerably smaller. The optical data for Coma cannot accept these lower 
values however, and the intersection of the two data sets suggest M/L « 100 to 
125 hgQ--'- for Coma. 

Typically, the virial mass of the clusters at 2 hgg - * Mpc lies around 
10 15 M q . For Coma, extrapolating the gas profiles to this radius would give a 
gas mass of around 2.10* 4 hgg" 5 / 2 M Q . The resldua1 mass 1n galaxies would 
correspond to a mass-to-light ratio of around 70 h 50 ~* if the gas were 
excluded. This remarkably low value would suggest that at least in the 
clusters, the dark mass of the galactic halos is not extremely large. 

Finally, the present result should encourage us to revisit our thinking 
on whether q = 1, since it removes one of the very few reasons to suppose that 
the mass-to-light ratio increases with scale size in the universe. 

We thank Keith Arnaud for help in analyzing the Coma velocity profile. 
R.M. thanks the Institute of Astronomy, Cambridge, England for hospitality. 
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FIGURE CAPTIONS 


1. The maximum virial mass as a function of radius for the five 
clusters. These values are obtained for the maximum delta core radius pair. 
These values are obtained from Jones and Forman 1984. 

2. a The distribution of temperature vs. radius for Coma for a delta = 1 
model. The solid line is for the case of M/Ly = 100 and the dotted line for 
M/L w = 150. 


2. b The distribution of temperature vs. radius for Coma for a delta = 1.25 
model. The solid line is for the case of M/Ly = 100 and the dotted line for 
M/Ly = 150. 

3. a The velocity dispersion vs. radius for Coma for the best fit polytropic 
models (solid line for the model that best fits the X-ray data, dotted line 
for the best fit to the optical data only) vs. the data points of Kent and 
Gunn. 


3.b The velocity dispersion vs. radius for the Perseus cluster for the best 
fit polytropic models (solid line for the model that best fits the X-ray data, 
dotted line for the best fit to the optical data only, vs. the data points of 
Kent and Sargent. 
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4. The mass to light ratio for Coma vs. radius using a delta = 1.0, gamma 
= 1.25, Tq = 15.3 keV model. Notice the rise in the center and the flatness 
of M/L vs. radius at large distances. The dashed line indicates a possible 
M/L value if the central 0.5 Mpc of Coma is better described by a power law 
(Tonry and Beers) galaxy distribution. 
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TEMPERATURE PROFILES 
(DELTA =1.0) 



Figure 2a 
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Figure 2b 
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Figure 4 





